Evolutionary origin and distribution of amino acid mutations associated with resistance to sodium channel modulators in onion thrips, Thrips tabaci.

In onion thrips Thrips tabaci, reduced sensitivity of the sodium channel caused by several sodium channel mutations have been correlated with pyrethroid resistance. For this study, using mitochondrial cytochrome c oxidase subunit I gene sequences, we examined the phylogenetic relation among a total of 52 thelytokous and arrhenotokous strains with different genotypes of the sodium channel mutations. Then, we used flow cytometry to estimate their ploidy. Results showed that the strains are divisible into three groups: diploid thelytoky, triploid thelytoky, and diploid arrhenotoky. Using 23 whole genome resequencing data obtained from 20 strains out of 52, we examined their genetic relation further using principal component analysis, admixture analysis, and a fixation index. Results showed that diploid and triploid thelytokous groups are further classifiable into two based on the sodium channel mutations harbored by the respective group members (strains). The greatest genetic divergence was observed between thelytokous and arrhenotokous groups with a pair of T929I and K1774N. Nevertheless, they shared a genomic region with virtually no polymorphism around the sodium channel gene loci, suggesting a hard selective sweep. Based on these findings, we discuss the evolutionary origin and distribution of the sodium channel mutations in T. tabaci.

A major endogenous glycoside hydrolase mediating quercetin uptake in Bombyx mori.

Quercetin is a common plant flavonoid which is involved in herbivore-plant interactions. Mulberry silkworms (domestic silkworm, Bombyx mori, and wild silkworm, Bombyx mandarina) take up quercetin from mulberry leaves and accumulate the metabolites in the cocoon, thereby improving its protective properties. Here we identified a glycoside hydrolase, named glycoside hydrolase family 1 group G 5 (GH1G5), which is expressed in the midgut and is involved in quercetin metabolism in the domestic silkworm. Our results suggest that this enzyme mediates quercetin uptake by deglycosylating the three primary quercetin glycosides present in mulberry leaf: rutin, quercetin-3-O-malonylglucoside, and quercetin-3-O-glucoside. Despite being located in an unstable genomic region that has undergone frequent structural changes in the evolution of Lepidoptera, GH1G5 has retained its hydrolytic activity, suggesting quercetin uptake has adaptive significance for mulberry silkworms. GH1G5 is also important in breeding: defective mutations which result in discoloration of the cocoon and increased silk yield are homozygously conserved in 27 of the 32 Japanese white-cocoon domestic silkworm strains and 12 of the 30 Chinese ones we investigated.

ABC transporter subfamily B1 as a susceptibility determinant of Bombyx mori larvae to Cry1Ba, Cry1Ia and Cry9Da toxins.

ATP binding cassette (ABC) transporters are a diverse family of transmembrane proteins. Specific subfamily members expressed in the lepidopteran midgut can act as susceptibility determinants for several insecticidal Bt Cry proteins. However, the susceptibility determinants to many Cry toxins still remain unclear. Therefore, we knocked out a series of ABC transporters that are highly expressed in the midgut of Bombyx mori larvae by transcription activator-like effector nuclease (TALEN)-mediated gene editing, and the lineages that became resistant to Cry toxins were searched by toxin overlay bioassay. As a result, the B. mori ABC transporter subfamily B1 (BmABCB1) knockout lineage showed 19.17-fold resistance to Cry1Ba, 876.2-fold resistance to Cry1Ia, and 29.1-fold resistance to Cry9Da, suggesting that BmABCB1 is the determinant of susceptibility to these toxins. BmABCC2 and BmABCC3 have been shown to be susceptibility determinants based on their function as receptors. Therefore, we next heterologously expressed these ABC transporters in HEK293T cells and performed a cell swelling assay to examine whether these molecules could exert receptor functions. As a result, BmABCB1-expressing cells showed swelling response to Cry1Ia and Cry9Da, and cells expressing PxABCB1, which is the Plutella xylostella ortholog of BmABCB1, showed swelling for Cry1Ba, suggesting that ABCB1 is a susceptibility determinant by functioning as a receptor to these toxins. Furthermore, in order to clarify how high binding affinity is based on receptor function, we performed surface plasmon resonance analysis and found that each KD of Cry1Ba, Cry1Ia, and Cry9Da to BmABCB1 were 7.69 × 10-8 M, 2.19 × 10-9 M, and 4.17 × 10-6 M respectively.

A male-killing gene encoded by a symbiotic virus of Drosophila.

In most eukaryotes, biparentally inherited nuclear genomes and maternally inherited cytoplasmic genomes have different evolutionary interests. Strongly female-biased sex ratios that are repeatedly observed in various arthropods often result from the male-specific lethality (male-killing) induced by maternally inherited symbiotic bacteria such as Spiroplasma and Wolbachia. However, despite some plausible case reports wherein viruses are raised as male-killers, it is not well understood how viruses, having much smaller genomes than bacteria, are capable of inducing male-killing. Here we show that a maternally inherited double-stranded RNA (dsRNA) virus belonging to the family Partitiviridae (designated DbMKPV1) induces male-killing in Drosophila. DbMKPV1 localizes in the cytoplasm and possesses only four genes, i.e., one gene in each of the four genomic segments (dsRNA1-dsRNA4), in contrast to ca. 1000 or more genes possessed by Spiroplasma or Wolbachia. We also show that a protein (designated PVMKp1; 330 amino acids in size), encoded by a gene on the dsRNA4 segment, is necessary and sufficient for inducing male-killing. Our results imply that male-killing genes can be easily acquired by symbiotic viruses through reassortment and that symbiotic viruses are hidden players in arthropod evolution. We anticipate that host-manipulating genes possessed by symbiotic viruses can be utilized for controlling arthropods.

Highly accurate genome assembly of an improved high-yielding silkworm strain, Nichi01.

The silkworm (Bombyx mori) is an important lepidopteran model insect and an industrial domestic animal traditionally used for silk production. Here, we report the genome assembly of an improved Japanese strain Nichi01, in which the cocoon yield is comparable to that of commercial silkworm strains. The integration of PacBio Sequel II long-read and ddRAD-seq-based high-density genetic linkage map achieved the highest quality genome assembly of silkworms to date; 22 of the 28 pseudomolecules contained telomeric repeats at both ends, and only four gaps were present in the assembly. A total of 452 Mbp of the assembly with an N50 of 16.614 Mbp covered 99.3% of the complete orthologs of the lepidopteran core genes. Although the genome sequence of Nichi01 and that of the previously reported low-yielding tropical strain p50T assured their accuracy in most regions, we corrected several regions, misassembled in p50T, in our assembly. A total of 18,397 proteins were predicted using over 95 Gb of mRNA-seq derived from 10 different organs, covering 96.9% of the complete orthologs of the lepidopteran core genes. The final assembly and annotation files are available in KAIKObase (https://kaikobase.dna.affrc.go.jp/index.html) along with a genome browser and BLAST searching service, which would facilitate further studies and the breeding of silkworms and other insects.

Reference Genome Sequences of the Oriental Armyworm, Mythimna separata (Lepidoptera: Noctuidae).

Lepidopteran insects are an important group of animals, including those used as biochemical and physiological model species in the insect and silk industries as well as others that are major agricultural pests. Therefore, the genome sequences of several lepidopteran insects have been reported. The oriental armyworm, Mythimna separata, is an agricultural pest commonly used to study insect immune reactions and interactions with parasitoid wasps as hosts. To improve our understanding of these research topics, reference genome sequences were constructed in the present study. Using long-read and short-read sequence data, de novo assembly and polishing were performed and haplotigs were purged. Subsequently, gene predictions and functional annotations were performed. To search for orthologs of the Toll and Immune Deficiency (IMD) pathways and for C-type lectins, annotation data analysis, BLASTp, and Hummer scans were performed. The M. separata genome is 682 Mbp; its contig N50 was 2.7 Mbp, with 21,970 genes and 24,452 coding sites predicted. All orthologs of the core components of the Toll and IMD pathways and 105 C-type lectins were identified. These results suggest that the genome data were of sufficient quality for use as reference genome data and could contribute to promoting M. separata and lepidopteran research at the molecular and genome levels.

Co-Expression Network and Time-Course Expression Analyses to Identify Silk Protein Regulatory Factors in Bombyx mori.

Bombyx mori is an important economic insect and an animal model in pharmacomedical research. Although its physiology has been studied for many years, the mechanism via which silk protein genes are regulated remains unclear. In this study, we performed two-step expression screening, namely co-expression network and time-course expression analyses to screen silk protein regulation factors. A co-expression network analysis using RNA-seq data that were obtained from various tissues, including the silk glands of B. mori, was performed to identify novel silk protein regulatory factors. Overall, 91 transcription factors, including some known ones, were found to be co-expressed with silk protein genes. Furthermore, time-course expression analysis during the fifth instar larvae stage revealed that the expression pattern of 13 novel transcription factors was highly relevant to that of silk protein genes and their known regulatory factor genes. In particular, the expression peak of several transcription factors (TFs) was detected before the expression of silk protein genes peak. These results indicated that a larger number of genes than expected may be involved in silk protein regulation in B. mori. Functional analyses of function-unknown transcription factors should enhance our understanding of this system.

A feasibility trial of genomics-based diagnosis detecting insecticide resistance of the diamondback moth.

BACKGROUND: Insecticide resistance management has been key for crop protection for over 70 years and is increasingly important because the development of new active ingredients has decreased in recent years. By monitoring the development of resistance in a timely manner, we can effectively prolong insecticide efficacy. Genomic-based diagnosis can reliably predict resistance development if information on resistant mutations against major pesticides is available. Here, we developed a feasibility trial of genomics-based diagnosis of insecticide resistance in diamondback moth (Plutella xylostella) populations in Nagano Prefecture, Japan. Amplicon sequencing analyses using a next-generation sequencer (Illumina MiSeq) for major insecticides, including diamides, pyrethroids, Bacillus thuringiensis (Bt) toxin (Cry1Ac), organophosphates, and spinosyns, were conducted. RESULTS: Mutations related to the resistance of pyrethroids, organophosphates, and diamides (flubendiamide and chlorantraniliprole) prevailed, while those of a diamide (cyantraniliprole), Bt (Cry1Ac), and spinosyns were scanty, suggesting that they are still effective. The results of the genomics-based diagnosis were generally concordant with the results of bioassays. Resistance development tendencies were generally uniform across Nagano. CONCLUSION: An insecticide-resistance management campaign can be conducted in Nagano Prefecture with a quick genomic-based diagnosis in early spring while bioassay is the only option for monitoring resistances whose mutations are unavailable. Our study is the first step in the future management of insecticide resistance in all significant pests. © 2022 Society of Chemical Industry.

Combination of target site mutation and associated CYPs confers high-level resistance to pyridaben in Tetranychus urticae.

Pyridaben is a mitochondrial electron transport complex I inhibitor. The H110R mutation in the PSST subunit has been reported as a major factor in pyridaben resistance in the two-spotted spider mite, Tetranychus urticae. However, backcross experiments revealed that the mutant PSST alone conferred only moderate resistance. In contrast, inhibition of cytochrome P450 (CYP) markedly reduces resistance levels in a number of highly resistant strains. It was reported previously that maternal factors contributed to the inheritance of pyridaben resistance in the egg stage, but the underlying mechanisms have yet to be elucidated. Here, we studied the combined effects of the PSST H110R mutation and candidate CYPs, as metabolic resistance factors, on pyridaben resistance in T. urticae. We found that the maternal effects of inheritance of resistance in the egg stage were associated with CYP activity. Analysis of differential gene expression by RNA-seq identified CYP392A3 as a candidate causal factor for the high resistance level. Congenic strains, where the alleles of both PSST and CYP392A3 were derived from a resistant strain (RR_i; i = 1 or 2) and a susceptible strain (SS_i) in a common susceptible genetic background, were constructed by marker-assisted backcrossing. RR_i showed upregulation of CYP392A3 and high resistance levels (LC50 > 10,000 mg L-1), while SS_i had LC50 < 10 mg L-1. To disentangle the individual effects of PSST and CYP392A3 alleles, we also attempted to uncouple these genes in RR_i. We conclude that given the variation in LC50 values and expression levels of CYP392A3 in the congenic and uncoupled strains, it is likely that the high pyridaben resistance levels are due to a synergistic or cumulative effect of the combination of mutant PSST and associated CYPs, including CYP392A3, but other yet to be discovered factors cannot be excluded.

Reference Transcriptome Data in Silkworm Bombyx mori.

Herein, we performed RNA-seq analysis of ten major tissues/subparts of silkworm larvae. The sequences were mapped onto the reference genome assembly and the reference transcriptome data were successfully constructed. The reference data provided a nearly complete sequence for sericin-1, a major silk gene with a complex structure. We also markedly improved the gene model for other genes. The transcriptomic expression was investigated in each tissue and a number of transcripts were identified that were exclusively expressed in tissues such as the testis. Transcripts strongly expressed in the midgut formed tight genomic clusters, suggesting that they originated from tandem gene duplication. Transcriptional factor genes expressed in specific tissues or the silk gland subparts were also identified. We successfully constructed reference transcriptome data in the silkworm and found that a number of transcripts showed unique expression profiles. These results will facilitate basic studies on the silkworm and accelerate its applications, which will contribute to further advances in lepidopteran and entomological research as well as the practical use of these insects.

Characterization of E93 in neometabolous thrips Frankliniella occidentalis and Haplothrips brevitubus.

Insect metamorphosis into an adult occurs after the juvenile hormone (JH) titer decreases at the end of the juvenile stage. This generally coincides with decreased transcript levels of JH-response transcription factors Krüppel homolog 1 (Kr-h1) and broad (br), and increased transcript levels of the adult specifier E93. Thrips (Thysanoptera) develop through inactive and non-feeding stages referred to as "propupa" and "pupa", and this type of distinctive metamorphosis is called neometaboly. To understand the mechanisms of hormonal regulation in thrips metamorphosis, we previously analyzed the transcript levels of Kr-h1 and br in two thrips species, Frankliniella occidentalis (Thripidae) and Haplothrips brevitubus (Phlaeothripidae). In both species, the transcript levels of Kr-h1 and br decreased in the "propupal" and "pupal" stages, and their transcription was upregulated by exogenous JH mimic treatment. Here we analyzed the developmental profiles of E93 in these two thrips species. Quantitative RT-PCR revealed that E93 expression started to increase at the end of the larval stage in F. occidentalis and in the "propupal" stage of H. brevitubus, as Kr-h1 and br mRNA levels decreased. Treatment with an exogenous JH mimic at the onset of metamorphosis prevented pupal-adult transition and caused repression of E93. These results indicated that E93 is involved in adult differentiation after JH titer decreases at the end of the larval stage of thrips. By comparing the expression profiles of Kr-h1, br, and E93 among insect species, we propose that the "propupal" and "pupal" stages of thrips have some similarities with the holometabolous prepupal and pupal stages, respectively.

Ovicidal activity of juvenile hormone mimics in the bean bug, Riptortus pedestris.

Insect juvenile hormone (JH) mimics (JHMs) are known to have ovicidal effects if applied to adult females or eggs. Here, we examined the effects of exogenous JHMs on embryonic development of the bean bug, Riptortus pedestris. The expression profiles of JH early response genes and JH biosynthetic enzymes indicated that JH titer was low for the first 3 days of the egg stage and increased thereafter. Application of JH III skipped bisepoxide (JHSB3) or JHM on Day 0 eggs when JH titer was low caused reduced hatchability, and the embryos mainly arrested in mid- or late embryonic stage. Application of JHMs on Day 5 eggs also resulted in an arrest, but this was less effective compared with Day 0 treatment. Interestingly, ovicidal activity of synthetic JHMs was much lower than that of JHSB3. This study will contribute to developing novel insecticides that are selective among insect species.

An update of KAIKObase, the silkworm genome database.

KAIKObase was established in 2009 as the genome database of the domesticated silkworm Bombyx mori. It provides several gene sets and genetic maps as well as genome annotation obtained from the sequencing project of the International Silkworm Genome Consortium in 2008. KAIKObase has been used widely for silkworm and insect studies even though there are some erroneous predicted genes due to misassembly and gaps in the genome. In 2019, we released a new silkworm genome assembly, showing improvements in gap closure and covering more and longer gene models. Therefore, there is a need to include new genome and new gene models to KAIKObase. In this article, we present the updated contents of KAIKObase and the methods to generate, integrate and analyze the data sets. Database URL: https://kaikobase.dna.affrc.go.jp.

Origin, selection, and spread of diamide insecticide resistance allele in field populations of diamondback moth in east and southeast Asia.

BACKGROUND: The investigation of molecular mechanisms and evolution of resistance to insecticides is an ongoing challenge, as researchers must provide guidance to manage the resistance to achieve sustainable production in agriculture. Predicting, monitoring, and managing insecticide resistance requires information on the origins, selection, and spread of resistance genes. The resistance of Plutella xylostella (L.) against diamide insecticides is becoming an increasingly severe problem in east and southeast Asia. In this study, the evolution of resistance was investigated using a resistance allele [ryanodine receptor (RyR); G4946E mutation] and its flanking regions, as well as mitochondrial cytochrome c oxidase subunit I (mtCOI). RESULTS: The sequences of the flanking region of the G4946E and mtCOI suggested that the G4946E mutation has a key role in resistance. Furthermore, the G4946E mutation has multiple origins, and congenic resistant mutations have spread across east and southeast Asia, despite substantial geographical barriers. In addition, the susceptibility of field populations partially recovered during winter, based on the observed decrease in the G4946E (resistant allele) frequency. Finally, the resistance level indexed by the frequency of the E4946 allele was significantly lower in non-overwintering regions than in overwintering regions. CONCLUSION: The information of the present study is useful to monitor resistance using molecular markers and to develop strategies to delay the evolution of diamide resistance.

Effectiveness of orally-delivered double-stranded RNA on gene silencing in the stinkbug Plautia stali.

Development of a reliable method for RNA interference (RNAi) by orally-delivered double-stranded RNA (dsRNA) is potentially promising for crop protection. Considering that RNAi efficiency considerably varies among different insect species, it is important to seek for the practical conditions under which dsRNA-mediated RNAi effectively works against each pest insect. Here we investigated RNAi efficiency in the brown-winged green stinkbug Plautia stali, which is notorious for infesting various fruits and crop plants. Microinjection of dsRNA into P. stali revealed high RNAi efficiency-injection of only 30 ng dsRNA into last-instar nymphs was sufficient to knockdown target genes as manifested by their phenotypes, and injection of 300 ng dsRNA suppressed the gene expression levels by 80% to 99.9%. Knockdown experiments by dsRNA injection showed that multicopper oxidase 2 (MCO2), vacuolar ATPase (vATPase), inhibitor of apoptosis (IAP), and vacuolar-sorting protein Snf7 are essential for survival of P. stali, as has been demonstrated in other insects. By contrast, P. stali exhibited very low RNAi efficiency when dsRNA was orally administered. When 1000 ng/µL of dsRNA solution was orally provided to first-instar nymphs, no obvious phenotypes were observed. Consistent with this, RT-qPCR showed that the gene expression levels were not affected. A higher concentration of dsRNA (5000 ng/µL) induced mortality in some cohorts, and the gene expression levels were reduced to nearly 50%. Simultaneous oral administration of dsRNA against potential RNAi blocker genes did not improve the RNAi efficiency of the target genes. In conclusion, P. stali shows high sensitivity to RNAi with injected dsRNA but, unlike the allied pest stinkbugs Halyomorpha halys and Nezara viridula, very low sensitivity to RNAi with orally-delivered dsRNA, which highlights the varied sensitivity to RNAi across different species and limits the applicability of the molecular tool for controlling this specific insect pest.

Transcriptomic analysis of the bagworm moth silk gland reveals a number of silk genes conserved within Lepidoptera.

Lepidopteran insects produce cocoons with unique properties. The cocoons are made of silk produced in the larval tissue silk gland and our understanding of the silk genes is still very limited. Here, we investigated silk genes in the bagworm moth Eumeta variegata, a species that has recently been found to produce extraordinarily strong and tough silk. Using short-read transcriptomic analysis, we identified a partial sequence of the fibroin heavy chain gene and its product was found to have a C-terminal structure that is conserved within nonsaturniid species. This is in accordance with the presence of fibroin light chain/fibrohexamerin genes and it is suggested that the bagworm moth is producing silk composed of fibroin ternary complex. This indicates that the fibroin structure has been evolutionarily conserved longer than previously thought. Other than fibroins we identified candidates for sericin genes, expressed strongly in the middle region of the silk gland and encoding serine-rich proteins, and other silk genes, that are structurally conserved with other lepidopteran homologues. The bagworm moth is thus considered to be producing conventional lepidopteran type of silk. We further found a number of genes expressed in a specific region of the silk gland and some genes showed conserved expression with Bombyx mori counterparts. This is the first study allowing comprehensive silk gene identification and expression analysis in the lepidopteran Psychidae family and should contribute to the understanding of silk gene evolution as well as to the development of novel types of silk.

Genome-wide assessment and development of molecular diagnostic methods for imidacloprid-resistance in the brown planthopper, Nilaparvata lugens (Hemiptera; Delphacidae).

BACKGROUND: The brown planthopper, Nilaparvata lugens (Stål), is one of the most notorious pests of rice throughout Asia. The brown planthopper has developed high resistance to imidacloprid, a member of neonicotinoid insecticides. Several genes and mutations conferring imidacloprid resistance in N. lugens, especially in eastern and southeastern Asia populations, have been reported. Thus, the key mechanisms of imidacloprid resistance need to be examined. RESULTS: RNA-seq analyses revealed that only one cytochrome P450 monooxygenase gene, CYP6ER1, was commonly upregulated in the five resistant strains tested. Sequences of CYP6ER1, which were highly expressed in the imidacloprid-resistant strains, contained a three-nucleotide deletion in the coding region, and amino acid substitutions and deletion, compared to that in an imidacloprid-susceptible strain. RNAi-mediated gene knockdown of CYP6ER1 increased imidacloprid susceptibility in a resistant strain. Further, we established two simple and convenient PCR-based molecular diagnostic methods to detect the CYP6ER1 locus with the three-nucleotide deletion. Using these methods, the resistance of F2 progenies derived from the crosses of F1 siblings from susceptible and resistant parents was analyzed, showing that the imidacloprid resistance had a relationship to the CYP6ER1 locus with the three-nucleotide deletion. CONCLUSION: The overexpression of a variant CYP6ER1 with amino acid substitutions and deletion was involved in imidacloprid resistance in N. lugens. Based on these findings, molecular diagnostic methods have been developed and are promising tools for monitoring imidacloprid resistance in paddy fields. © 2020 Society of Chemical Industry.

QTL mapping using microsatellite linkage reveals target-site mutations associated with high levels of resistance against three mitochondrial complex II inhibitors in Tetranychus urticae.

The acaricides cyflumetofen, cyenopyrafen, and pyflubumide act as inhibitors of the mitochondrial electron transport system at complex II (succinate dehydrogenase; SDH), a new mode of action in arthropods. The development and mechanisms of low-level resistance against cyenopyrafen and cyflumetofen have been previously reported in Tetranychus urticae. In the present study, we investigated high levels of resistance against three SDH inhibitors in T. urticae field populations and clarify the genetic basis of resistance using quantitative trait locus (QTL) analysis. First, we constructed a microsatellite linkage map comprising 64 markers assembled into three linkage groups (LGs) with total length of 683.8 cM and average marker spacing of 11.03 cM. We then used the linkage map to perform QTL mapping, and identified significant QTLs contributing to resistance to cyflumetofen (one QTL on LG1), cyenopyrafen (one QTL on LG3), and pyflubumide (two QTLs on LG1 and LG3). The QTL peaks on LG1 for cyflumetofen and pyflubumide overlapped and included the SdhB locus. For cyenopyrafen resistance, the QTLs on LG3 included the SdhC locus. For cyflumetofen resistance, we found an I260T mutation in SdhB. For pyflubumide and cyenopyrafen resistance, we detected I260V and S56L substitutions in SdhB and SdhC, respectively, by direct sequencing. Both I260 in SdhB and S56 in SdhC were present in highly conserved regions of the ubiquinone binding site formed at the interface among SdhB, SdhC, and SdhD. Mutations at these positions have been implicated in resistance against fungicides that act as Sdh inhibitors in various pathogens. Therefore, we consider these mutations to be target-site resistance mutations for these acaricidal SDH inhibitors.

Sex-specific expression profiles of ecdysteroid biosynthesis and ecdysone response genes in extreme sexual dimorphism of the mealybug Planococcus kraunhiae (Kuwana).

Insect molting hormone (ecdysteroids) and juvenile hormone regulate molting and metamorphic events in a variety of insect species. Mealybugs undergo sexually dimorphic metamorphosis: males develop into winged adults through non-feeding, pupa-like stages called prepupa and pupa, while females emerge as neotenic wingless adults. We previously demonstrated, in the Japanese mealybug Planococcus kraunhiae (Kuwana), that the juvenile hormone titer is higher in males than in females at the end of the juvenile stage, which suggests that juvenile hormone may regulate male-specific adult morphogenesis. Here, we examined the involvement of ecdysteroids in sexually dimorphic metamorphosis. To estimate ecdysteroid titers, quantitative RT-PCR analyses of four Halloween genes encoding for cytochrome P450 monooxygenases in ecdysteroid biosynthesis, i.e., spook, disembodied, shadow and shade, were performed. Overall, their expression levels peaked before each nymphal molt. Transcript levels of spook, disembodied and shadow, genes that catalyze the steps in ecdysteroid biosynthesis in the prothoracic gland, were higher in males from the middle of the second nymphal instar to adult emergence. In contrast, the expression of shade, which was reported to be involved in the conversion of ecdysone into 20-hydroxyecdysone in peripheral tissues, was similar between males and females. These results suggest that ecdysteroid biosynthesis in the prothoracic gland is more active in males than in females, although the final conversion into 20-hydroxyecdysone occurs at similar levels in both sexes. Moreover, expression profiles of ecdysone response genes, ecdysone receptor and ecdysone-induced protein 75B, were also analyzed. Based on these expression profiles, we propose that the changes in ecdysteroid titer differ between males and females, and that high ecdysteroid titer is essential for directing male adult development.

Diversity and function of multicopper oxidase genes in the stinkbug Plautia stali.

Multicopper oxidase (MCO) genes comprise multigene families in bacteria, fungi, plants and animals. Two families of MCO genes, MCO1 (laccase1) and MCO2 (laccase2), are conserved among diverse insects and relatively well-characterized, whereas additional MCO genes, whose biological functions have been poorly understood, are also found in some insects. Previous studies reported that MCO1 participates in gut immunity and MCO2 plays important roles in cuticle sclerotization and pigmentation of insects. In mosquitoes, MCO2 was reported to be involved in eggshell sclerotization and pigmentation, on the ground that knockdown of MCO2 caused deformity and fragility of the eggshell. Here we identified a total of 7 MCO genes, including PsMCO1 and PsMCO2, and investigated their expression and function in the brown-winged green stinkbug Plautia stali. RNA interference (RNAi) knockdown of MCO genes by injecting double-stranded RNA (dsRNA) into nymphs revealed that MCO2, but not the other 6 MCOs, is required for cuticle sclerotization and pigmentation, and also for survival of P. stali. Trans-generational knockdown of MCO2 by injecting dsRNA into adult females (maternal RNAi) resulted in the production of unhatched eggs despite the absence of deformity or fragility of the eggshell. These results suggested that MCO2 plays an important role in sclerotization and pigmentation of the cuticle but not in eggshell integrity in P. stali. Maternal RNAi of any of the other 6 MCO genes and 3 tyrosinase genes affected neither survival nor eggshell integrity of P. stali. Contrary to the observations in the red flour beetle and the brown rice planthopper, RNAi knockdown of MCO6 (MCORP; Multicopper oxidase related protein) exhibited no lethal effects on P. stali. Taken together, our findings provide insight into the functional diversity and commonality of MCOs across hemipteran and other insect groups.